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The Sunyaev-Zel'dovich (SZ) effect expected to be associated with massive star- 
formation activity produced during the formation of the most luminous bulges of 
normal galaxies is discussed. 

Using 1-D chemohydrodynamical models for spheroidal galaxy evolution we show 
that during the early epochs of galaxy evolution the gas in massive events of star 
formation may reach temperatures and densities high enough to produce values of the 
comptonisation parameter y comparable to those present in galaxy clusters. 

In this scenario, we discuss the possibility of detection of the SZ signature in high 
redshift starforming galaxies with the next generation of mm telescopes capable of 
arcsecond resolution and equipped with high sensitivity detectors. 

We show how millimeter colour-colour diagrams or diagnostic diagrams could be 
used to distinguish between the dust emission and the SZ effect and suggest the use 
of simultaneous multifrequency observations to improve the chances of detecting the 
SZ effect. 
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1 INTRODUCTION 

The SZ effect (Sunyaev & Zel'dovich 1972) is produced when 
the cosmic microwave background (CMB) photons from the 
Rayleigh- Jeans region, due to the interaction with the fast 
electrons, move to the Wien tail of the CMB spectrum thus 
producing a unique spectral signature. This signature shows 
an increase with respect to the mean CMB brightness of the 
observed intensity for wavelengths shorter than 1.34 mm 
and a decrease for wavelengths longer than 1.34 mm. This 
is known as the thermal SZ effect. Another component is 
the kinetic SZ effect, due to the interaction with CMB pho- 
tons of plasma coherently moving with respect to the refer- 
ence frame where the CMB is isotropic (e.g. Rephaeli 1995, 
Birkinshaw 1999, Church, Jaffe and Knox 2001). 

The SZ effect constitutes a powerful tool to study the 
physical properties of the intercluster gas, and when com- 
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bined with X-ray observations, also to estimate cosmological 
constants such as the Hubble constant Ho, the current mass 
density of the universe Qm or the cosmological constant A 
(e.g. Majumdar 2001, Diego et al. 2002, Levine, Schulz & 
White 2002, Reese et al. 2002, Majumdar & Mohr 2003). 

It is understood that in clusters the thermal SZ effect 
due to the interaction of CMB photons with electrons with 
temperatures around 10 8 K dominates over the kinetic SZ 
effect. This fact is evidenced by combining observations at 
1 mm where the thermal effect is close to the maximum, at 
1.4 mm where it is almost zero and the kinetic SZ effect can 
be marginally detected and at 2 mm where the thermal SZ 
has its minimum and is negative (e.g. Lamarre et al. 1998, 
Mauskopf et al. 2000). 

The advent of a new generation of mm telescopes 
(LMT/GTM or ALMA, described in section 3) that com- 
bine relatively high angular resolution with high sensitivity 
will allow to explore the possibility of detecting the SZ effect 
in volumes much smaller than that of a cluster of galaxies. 
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Would it be possible for these new instruments to detect the 
SZ signature in individual galaxies? 

Natarajan &: Sigurdsson (1999), Rosa-Gonzalez et 
al. (2000, 2001), Majumdar, Nath and Chiba (2001) and 
Rosa-Gonzalez (2002) have indeed explored the detactabil- 
ity of the SZ effect in individual galaxies. 

Natarajan & Sigurdsson (1999) studied the SZ effect 
due to galactic outflows powered by the mechanical energy 
provided by the accretion of matter to the central supermas- 
sive black hole in QSOs. Majumdar et al. (2001) propose 
the evolution of supernova driven galactic winds during the 
early stages of evolution of a normal galaxy, as the cause of 
the distortion of the CMB radiation due to the kinetic SZ 
effect. A similar model was originally proposed by Tegmark, 
Silk and Evrard (1993) to explain the absence of absorp- 
tion lines in the observed spectra of high redshift quasars 
(Gunn-Peterson effect, Gunn & Peterson 1965). The explo- 
sion model assumes an initial energy input which is equal 
to the 2% of the total luminosity generated in the super- 
nova explosions during about 5xl0 7 years. The range of 
galaxy masses covered by this model goes from 5xlO 7 M0 
to lO n M . 

Rosa-Gonzalez et al. (2000, 2001) and Rosa-Gonzalez 
(2002) discussed the possibility of observing the SZ effect 
in the centre of young giant spheroidal galaxies. They in- 
clude in the calculations both the thermal and the kinetic 
SZ effect. 

In the accepted galaxy formation scenario described by 
the hierarchical model (e.g. Navarro, Frenk & White 1995), 
massive galaxies form after merging with smaller system 
that have been formed in a previous epoch. However, there 
is observational evidence that a significant number of L* 
galaxies exist at redshift of about 5 and beyond (Hu, Cowie 
& McMahon 1998, Stanway, Bunker & McMahon 2003). 
Jimenez et al. (1999) discussed the possibility of the exis- 
tence of massive elliptical galaxies at high redshift, contrary 
to early predictions of the hierarchical models of galaxy for- 
mation. The study of massive galaxies at high redshift that 
could be detected through the unique SZ signature could be 
used to constrain the different galaxy formation scenarios. 

The SZ effect due to young galaxies together with radio 
sources (Holder 2002) and dusty galaxies (Blain, Ivison & 
Smail 1998, Blain 1998) are potential sources of noise in sen- 
sitive studies of SZ clusters and upcoming measurements of 
the CMB fluctuations (e.g. FIRST, Plank). Therefore, their 
characteristic emission and physical scale must be quanti- 
fied. 

In this paper we explore in detail the possibility of de- 
tection of the SZ effect in young galaxies undergoing a mas- 
sive starformation event. We discuss the confusion due to 
dust emission and suggest a diagnostic to segregate between 
dust emission and SZ effect distortions. 
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where the energy of the CMB photons is given by hv , the 
energy of the photons which come out from the cluster is 
hv, s is defined by s = log(u/uo), r e is the optical depth 
for scattering and P(s) is the probability that an incident 
photon with energy hv becomes a photon with an energy 
hv (Birkinshaw 1999). 

The integral in the Equation 1 can be calculated by us- 
ing the non relativistic expression given by the Kompaneets 
approximation (Kompaneets 1957), 
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where the distortion of the photon spectrum, h(x) is defined 
by 

h ( x ) = x (e /_ 1)2 ( x coth f - 4 ) ( 3 ) 

x is the dimensionless frequency, x = hv/(kT rad ) with 
Trad = 2.728 the temperature of the CMB, 
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and the Comptonization parameter "y " is defined by the 
integral of the electron pressure along the line of sight 
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In this equation gt is the Thomson scattering cross section, 
n e and T e are the electron number density and the electron 
temperature respectively, ks is the Boltzmann constant and 
m e c 2 is the rest energy of an electron; (see Rybicki & Light- 
man 1980 for a more detailed derivation of the different ex- 
pressions described above and Birkinshaw 1999 for a recent 
review). The effect of the hot gas on the CMB spectrum is to 
produce a fractional change in the brightness of the CMB 
radiation. We remark that the amplitude of the SZ effect 
is redshift-independent and depends only on the integrated 
properties of all the hot gas along the line of sight. The SZ 
effect can be considered as a tool to detect high pressure 
electron gas over cosmological distances. 



1.2 Kinetic SZ effect 

The kinetic SZ effect, as already mentioned, is due to the 
movement of the plasma with respect to the CMB rest 
frame. For a shell of matter which is moving with a radial 
velocity v with respect to the Hubble flow, the change of the 
intensity due to the kinetic SZ effect (Alkin) is given by, 



1.1 Thermal SZ effect 

The SZ effect describes the brightness change of the CMB 
due to the inverse Compton scattering by hot electrons (Sun- 
yaev & Zel'dovich 1972). 

The change in intensity between the unaffected CMB 
radiation and the radiation which go through the hot gas 
can be expressed by the integral, 
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where r e is the scattering optical depth, lo is given by equa- 
tion 4 and x is the dimensionless frequency. 

At low frequencies, the ratio between the kinetic SZ 
effect (Alkin) to the thermal SZ effect (AIth) is independent 
of the scattering optical depth and can be calculated by, 

AI kin 1 v m e c 2 ^ 
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Figure 2. Electron density and temperature profiles for the model of baryonic mass 2xfO 12 M0. The different lines correspond to different 
times: 20 Myr for the solid line, 600 Myr for the dotted line, 1 Gyr for the dashed line, 4 Gyr for the dot-dashed line and f f Gyr for the 
dot-dot-dashed line. 
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Figure 1. The kinetic SZ effect (solid and dot-dashed lines) sig- 
nature is compared to the thermal SZ effect (dashed line) for the 
case where the kinetic effect is half of the thermal one. The two 
cases of the kinetic SZ effect correspond to gas moving away from 
(dot-dashed line) and towards (solid line) the observer. 

where v is the velocity of the flow, m e c 2 is the electron rest 
energy and T e is the gas temperature (Birkinshaw 1999). 
The signatures of the thermal and kinetic SZ effects are 
plotted in Figure 1. 



2 THE GALAXY FORMATION MODEL 

In order to estimate the amplitude and time evolution of the 
SZ effect for individual galaxies, it is necessary to know the 
time and radial dependence of the gas temperature, density 
and velocity. We have restricted our study to those stages 
of galaxy evolution characterized by a high injection rate 
of mass and energy into the galaxy ISM (which happens to 
occur during the first few Gyr of the galaxy evolution). We 
have further assumed that this occurs during a maximum 
in the star formation rate produced by any of the known 
mechanisms. 

The evolution and radial profiles mentioned above (Fig- 



ure 2) have been obtained from a chemohydrodynamical 
model for evolution of elliptical galaxies (Friaga & Ter- 
levich 1998, hereafter FT98). The model combines multi- 
zone chemical evolution with 1-D hydrodynamics to calcu- 
late the evolution of chemical abundances in gas and stars, 
and also the temperature, density and velocity profiles of 
the gas. In this way, we can calculate the y profiles which 
are used to estimate the flux collected by a given beam size 
at a given frequency. 

In the FT98 model, a single massive dark halo hosts 
baryonic gas that falls toward the centre of the dark halo 
and will subsequently form stars. The dark halo is given by a 
static mass density distribution ph{r) = Pho\\ + (r/r fc ) 2 ] -1 , 
where pho is the halo central density and tk is the halo 
core radius. The gas and the stars exchange mass through 
star formation and stellar mass losses (supernovae, plane- 
tary nebulae, and stellar winds). Both the stellar distribu- 
tion and the dark halo are truncated at a common tidal ra- 
dius r t . The system, assumed to be spherical, is subdivided 
in several spherical zones and the hydrodynamical evolution 
of its ISM is calculated. The equations of chemical evolu- 
tion for each zone are then solved taking into account the 
gas flow, and the evolution of the chemical abundances is 
obtained. A total of ^100 star generations are stored dur- 
ing 13 Gyr for chemical evolution calculations. We assume 
that at a given radius r and time t, the specific star forma- 
tion rate v{r,t) follows a power-law function of gas density 
(p): v(r,t) = v{p/ Po) 1 ^ 2 , where po is the initial average gas 
density inside ru, and v is the normalization of the star 
formation law. We include inhibition of star formation for 
expanding gas (V.u > 0) or when the density is too low, 
and, therefore, the cooling is inefficient (i.e., for a cooling 
time t coo = (3/2)fcsT//im_r/A(T)p longer than the dynami- 
cal time tdyn = (37r/16 G p) 1 2 ) by multiplying v as defined 
above by the inhibition factors (1 + td y n max(0, V.tt)) -1 and 
(1 + tcoo/tdyn)^ 1 . A characteristic of these models is that 
several episodes of inflow and outflow occur simultaneously 
at different radii. The chemodynamical model for spheroids 
was used to investigate the relation between young elliptical 
galaxies and QSO activity (FT98), the absence of passively 
evolving elliptical galaxies in deep surveys (Jimenez et al. 
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1999), Lyman Break Galaxies (Friaca & Terlevich 1999), 
Blue Core Spheroids (Friaga & Terlevich 2001), the coupled 
spheroid and black hole formation (Archibald et al. 2002), 
and the link between DLAs and dwarf galaxies (Lanfranchi 
& Friaca 2003). 

The models are parameterized according to the (initial) 
baryonic mass inside the tidal radius, Mg = M g + M*, rn, 
r t , and the ratio of the halo to the (initial) luminous mass, 
Mh/Mq- We have investigated a grid of runs with Mg be- 
tween 10 11 and 5 x 10 12 Mq, and rn in the range 2.5 — 15 
kpc (see FT98). We set r t = 28r h and M h /M G = 5. That 
value is compatible with the recent WMAP results where the 
baryon density f2b = 0.047 ± 0.006 and the matter density 
Q. M = 0.29 ± 0.07 (Spergel et al. 2003). Smaller (larger) 
values of the ratio dark matter versus baryonic mass do not 
significantly alter the observed fluxes, but due to the shal- 
lower (deeper) potential well, the galactic winds will appear 
earlier (later). For example, in the model with Mq = 10 12 
Mq, the galactic wind occurs at 1.51, 1.55 and 1.67 Gyr, for 
M h /MG = 3, 5 and 7.5 respectively. 

The chemical evolution is driven by the stellar winds, 
planetary nebulae and SN phase which produce the enrich- 
ment of the ISM. New generation of stars form in the more 
metal-rich medium. We do not assume instantaneous recy- 
cling approximation for the chemical enrichment, but in- 
stead we take into account the delays for gas restoring from 
the stars due to the main-sequence lifetimes. Instantaneous 
mixing with the ISM is assumed for the stellar ejecta. We 
use metallicity dependent yields for SNe II, SNe la, and in- 
termediate mass stars (IMS): the SNe II yields of Woosley & 
Weaver (1995), for metallicities Z/Z = 0, 10 -4 , 10~ 2 , 10 _1 
and 1; SNIa yields from Iwamoto et al. (1999) — their mod- 
els W7 and W70 (progenitor metallicity Z=Zq and Z=0, 
respectively); the yields for IMS (0.8 — 8 M ), with initial 
Z=0.001, 0.004, 0.008, 0.02 and 0.4, from van den Hoek & 
Groenewegen (1997) (their variable t]agb case). For more 
details of the nucleosynthesis prescriptions, see FT98. The 
models start with an entirely gaseous protogalaxy with pri- 
mordial chemical abundances (Y — 0.24, Z = 0). 

In the FT98 models there is self-consistency of the hy- 
drodynamics, chemical evolution and atomic physics, since 
the cooling function is evaluated based on the actual chemi- 
cal abundances obtained from the chemo-dynamical mod- 
elling. For the sake of simplicity, our cooling function is 
a function only of the abundances of O and Fe, which 
are the main coolants for T > 10 J K. In the calculation 
of the cooling function, the abundances of elements other 
than Fe and O have been scaled to the O abundance as 
rii — rii t p + (n^Q — rii t p)no/no,Q, where rii is the ISM 
abundance by number of the element i, rii } p its primordial 
abundance and n^,© its solar abundance [the photospheric 
values of Holweger (2001) for N and O , and the meteoritic 
values of Grevesse & Sauval (1998) for the other elements]. 

It is important to note that the FT98 models describe 
the evolution of isolated galaxies and do not take into ac- 
count effects such as mergers that play an important role 
in structure formation within hierarchical galaxy formation 
models. However, the FT98 models can be considered as a 
description of the evolution of the hot gas after a major star 
forming event produced by interaction or merger. In addi- 
tion, while the present chemodynamical model accounts for 
complex flow patterns, with inflow in some regions and out- 



flow in other regions, the model does not take into account 
asymmetric fluxes due to the assumed spherical symmetry. 
Despite the limitations due to the spherical symmetry, the 
detailed predictions of this model are useful due to the good 
radial resolution of the galaxy and the realistic treatment 
of the star formation and proper calculation of the chemical 
evolution of the gas and stars. The FT98 model represents 
a complementary tool to the semi-analytic models of galaxy 
formation, which currently adopt very simplified recipes of 
star-formation but with coarser space resolution (e.g. Som- 
merville & Primack 1999, Benson et al. 2001). 

An important aspect relevant to the present calculation 
is that the models provide an estimate of the size of the 
hot gas region. This is central for deriving the expected flux 
because the effective y parameter (which is the average of the 
y parameters within a given beam size, Section 2.1) depends 
strongly on the size of the telescope beam. 

2.1 The SZ effect in massive young galaxies. 

The FT98 models predict that the hot gas in massive star 
forming regions with ages between 10 s and 10 9 years can 
reach temperatures exceeding 10 7 K and densities higher 
than 1 cm~ 3 within a radius of a few kpc. These tempera- 
tures and densities imply values of the Comptonization pa- 
rameter y larger than 10~ 4 . The maximum temperatures 
of the gas in the models (~ 1 keV) are well below the 
relativistic limit (~ 500 keV), and, therefore, we use the 
non-relativistic approximation (e.g. Birkinshaw 1999) to es- 
timate the SZ effect (Equation 2). Notice that the central y 
parameter is calculated for a region of a few hundred parsecs 
across. Only in this region is the y parameter comparable to 
those observed in galaxy clusters. 

Figure 3 shows the variation of the thermal parameter 
(yth) with radius. In the beginning of the calculations (the 
t = 20 Myr lines in Figures 2 and 3) , the gas has not yet set- 
tled in the potential well of the galaxy and its density is low 
and temperature is cold everywhere (the model initial con- 
ditions assume gas at T = 10 4 K), and, as a consequence, y 
is very small. As the model evolves, the gas falls towards the 
centre and is compressed, giving rise to shocks that rapidly 
heat the gas in the core to approximately the virial tempera- 
ture of the system (T ~ 10 7 K). Then, a highly efficient star 
formation is occurring throughout the galaxy and a young 
stellar population is rapidly built up. Following the initial 
violent star formation burst, the first Type II SNe appear 
and heat the ISM. The star formation is vigorous and in- 
tense during the first ~ 0.5 Gyr, when nearly half of the 
stellar population is formed. In Figures 2 and 3, we see that 
at t = 0.6 Gyr, there is enough gas at temperatures around 
T = 10 7 K to maintain values of y close to 10~ 4 in the cen- 
tral kpc. At 1 Gyr since the start of the calculations, the 
star formation gas consumed most of the gas and the SNII 
(and SNIa) heating is past its maximum. Although the gas 
temperature is higher than before, its density has dropped 
along most of the galaxy, resulting in a significant decrease 
of y . As the galaxy evolves, the star formation rate de- 
creases, the main heating source of the ISM is provided by 
Type la supernovae. Eventually, the thermal energy of the 
gas is enough to overcome the escape velocity, and the re- 
maining gas is expelled from the galaxy by galactic winds. 
At this point, the primary burst of star formation ceases. 
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Figure 3. The thermal y parameter as a function of radius. Dif- 
ferent lines correspond to different times: 20 Myr for the solid 
line, 600 Myr for the dotted line, 1 Gyr for the dashed line, 4 Gyr 
for the dot— dashed line and 11 Gyr for the dot-dot-dashed line. 
The baryonic mass of the galaxy is 2x10 12 Mq. 
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Figure 4. Angular size corresponding to a physical size of 1 kpc 
for different cosmologies. The solid line is the result for a universe 
with Qm = 0.3 and f2yv = 0.7. The dotted line is for the case of 
S!m = 1 and Qa = and the dashed line is for Qm = 0.05 and 
U A = 0. 

The curves at t = 4 and 11 Gyr illustrate the post-galactic 
wind stage, when the gas temperature is high, but the low 
density implies a low y. 

The y parameter is smaller in the external layers where 
the pressure is lower, therefore the resulting y parameter 
averaged within the telescope beam is expected to be lower 
than the central one. This effect can be quantified by an 
effective y parameter (y e ff) defined as the average of the 
y parameter in the different shells weighted by the corre- 
sponding projected area and convolved with the telescope 
beam profile, 

y*ff= (J*f yds ) * B{9) (8) 

where 9 is the size of the telescope beam and B(9) is the 
telescope beam profile. 

In reality, the telescope beam will cover a variable frac- 



tion of the galaxy depending on its intrinsic size, its redshift 
and the adopted cosmology (Figure 4). In order to translate 
physical sizes to angular sizes, we adopted a cosmological 
model given by Qm = 0.3, Qa = 0.7 and h=0.7. 

To illustrate the effect of beam dilution we show in Fig- 
ure 5 the change of y e ff with respect to the beam angle 
for models at 400 Myr. The results are similar for differ- 
ent times. The behaviour of y e ff can be analysed in three 
separate cases, according to the beam size: 

• Small beam, y e ff increases as 9 2 , where 9 is the size 
of the beam. In this case the beam size is smaller than the 
first radius of the model at the given time, so the value of y 
increases until the telescope beam reaches the inner radius 
of the model galaxy. 

• Medium beam, for beam angles similar to the inner 
radius, the model obtains the maximum value of y which 
corresponds to the central one (y c ). The integrated pressure 
in the external layer is lower than in the inner ones, hence 
increasing the beam size produces a slight decrease of y e ff- 

• Large beam, when the size of the beam is larger than 
the region where the values of y are significant the averaged 
y goes down as 9~ 2 . Notice that for the most massive galaxy 
(baryonic mass 5xl0 12 M Q ) the fast decline is reached after 
an angle of about 10" but for the low mass galaxies the fast 
decline starts when the angle is about 1". 

Figure 6 shows that for the massive galaxies, the comp- 
tonisation parameter reaches values higher than 3x10 . In 
the case of quasar outflows, Natarajan & Sigurdsson (1999) 
estimated a value for the comptonisation parameter of about 
5xl0~ 5 . Similar values of the comptonisation parameter 
were estimated for the case of early galactic winds (Mar 
jumdar et al. 2001). 

However, the dilution effect produces a significant de- 
crease of the y parameter and consequently a smaller ex- 
pected flux (Figure 6). Using the output of the models, we 
obtained the effective y parameter as a function of time 
and beam angle. The behaviour of the central y and the 
effective y parameter for the case of a beam of 5" is plot- 
ted in Figure 6. The central y parameters have the max- 
ima between 150 xlO 6 years and 10 9 years, independently of 
the baryonic mass of the galaxy. The maximum values go 
from 2xl0~ 4 for galaxies with baryonic masses of 5xl0 12 
and 2xl0 12 M Q down to 3xl0 -5 for galaxies with baryonic 
masses of 1O 11 M0. The effective y within an angle of 5"is 
about 10 times lower than the central one (Figure 6). The 
given times in Figures 6, 7 and 9 correspond to galaxies 
which started to collapse at redshift 50. Within this initial 
conditions, the maximum SZ amplitude is reached at red- 
shifts between 12 and 5. 

2.2 The kinetic SZ effect in young galaxies 

There are two processes that could produce a significant 
kinetic SZ effect in young galaxies. One is the infall of the 
primordial gas to the center of the galaxy and the other 
is the outflow produced by SN explosions and winds from 
massive young stars. 

The calculations of the thermal SZ effect presented in 
the previous section, are for a spherical galaxy (the mod- 
els by FT98 on which the calculations are based, assume 
spherical symmetry). For the case of the kinetic SZ effect 



© 2003 RAS, MNRAS 000, 1-11 



6 D. Rosa- Gonzalez et al. 





Log Time (years) Log Time (years) 

Figure 6. Evolution of the y parameter. The left panel is the evolution of the central thermal y parameter (j/c). The right panel is the 
effective thermal y parameter (j/ e //) assuming a beam size of 5". The different lines correspond to different galaxy masses, coded as in 
Figure 5. 



of a spherical galaxy with isotropic winds, the total contri- 
bution to the kinetic y parameter becomes zero. In reality, 
the geometry of galactic winds tend to be more bi-conical 
than spherical. If there is asymmetry between the two cones, 
there will be a net kinetic SZ effect. Even within the context 
of our spherically symmetric calculations, we can estimate 
the kinetic SZ effect, in the case of extreme asymmetry, by 
considering the contribution of only one hemisphere of the 
galaxy. 

In this gross approximation, the kinetic effect is domi- 
nant at early phases of the galaxy evolution (less than 0.5 
Gyr) due to the infall of gas to the centre of the galaxy. Dur- 
ing this epoch the gas is relatively cold with temperatures 
of about 10 4 K favouring the kinetic effect over the thermal 
one. During this early phase the SZ kinetic effect could be 
several times higher than the thermal effect (Figure 7). In 
the central region of the galaxy the thermal effect domi- 
nates due to the high temperatures and densities and as a 
result the central yth parameter is higher than the kinetic 
component. 



Notice, however, that an assumption of extremely asym- 
metric infall is probably less justified than that of asym- 
metric wind. However, when the galactic wind is operating 
(t > 1 Gyr) due to its low density, the implied kinetic SZ 
effect would be much less than the thermal SZ effect. More- 
over at these late times, the galaxy is far from the time of 
maximum SZ effect (t ~ 0.5 Gyr). 



3 TELESCOPES AND DETECTORS 

There are a few future mm and sub-mm facilities that may 
be capable of detecting the SZ effect produced by individual 
galaxies. 

One is the Atacama Large Millimeter Array (ALMA), 
a project lead by institutions of Europe, Japan and North 
America. ALMA will be an array of 64 antennas of 12 m 
in diameter getting baselines up to 10 km. It is expected 
that ALMA will cover the wavelength range from 8 mm to 
350 fim (Blain 1999). It will be situated in the Atacama 
desert in Chile. 
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Figure 7. The total flux (solid line) due to the sum of the thermal 
(dotted line) and kinetic (dashed line) SZ effects for a galaxy of 
2x10 12 Mq (baryonic mass). We assumed that the kinetic effect 
has a positive sign so the corresponding flux is added directly to 
the flux due to the thermal SZ. The beam angle (5") corresponds 
to observations at 1 mm with a 50 meter telescope. 

Another one is the Large Millimeter Telescope 
(LMT/GTM), a binational project between UMASS (USA) 
and INAOE (Mexico), with a single dish antenna of 50 m di- 
ameter. It will be situated on the Cerro la Negra in Mexico. 
In the case of observations of the continuum, that is the 
case of the SZ effect, the LMT/GTM will be equipped with 
BOLOCAM, a sensitive wide field camera that will oper- 
ate between 3 and 1 mm. At 1 mm on the LMT/GTM, 
BOLOCAM will have a projected pixel size of 5", three 
times smaller than SCUBA on the JCMT observing at 
850 iim (Holland et al. 1998). 

We have also included the characteristics of the Green 
Bank Telescope (GBT) a 100 m single dish which is planned 
to be able to work at 3 mm in the next few years (Dicker 
2002). 

The mapping speed is defined as the ratio between the 
field of view (FOV) and the square of the sensitivity. High 
mapping speed is necessary to complete surveys covering big 
portions of the sky and to detect the fainter sources. 

Table 1 shows the sensitivity, FOV and mapping speed 
for the different telescopes. The sensitivity of BOLOCAM 
in the LMT/GTM is estimated from Glenn et al. 1998. The 
sensitivity and FOV of ALMA are from Blain (1999). Notice 
that ALMA is not designed to perform wide area surveys (so 
it does not have an associated mapping speed). 

We are interested in the detection of the SZ effect in 
young massive galaxies that at high redshift may have a 
few arcseconds in size and an estimated flux of several /xjy 
(See Section 2.1). Table 1 shows that only with the next 
generation of mm telescopes it may be possible to detect 
them. 



4 THE EFFECT OF DUST EMISSION 

Dust emission is a dominant component in the observed 
multifrequency spectrum of galaxies over a wide range of 
redshifts. Given the small signal associated with the SZ of 
individual galaxies, dust emission will probably dominate 



the mm spectrum making difficult the detection of the SZ 
signature. At high redshift the peak of the dust emission 
(~100 fim) moves to sub-mm and mm wavelengths (e.g. 
Rowan-Robinson et al. 1997, Hughes et al. 1998) i.e. over 
the region where the SZ effect has its maximum. 

We estimated the relative amplitude of the flux due to 
the SZ effect and that due to the dust by comparing the mm 
flux due to dust emission in M82, a well known starforming 
galaxy, with the expected flux due to the SZ effect. 

M82 is a nearby (~ 3.2 Mpc) starburst galaxy with an 
infrared luminosity of L(8 — lOOO^im) = 3 x 10 10 L Q which 
represents the 84% of the total galaxy luminosity (Telesco 
1988). Hughes, Gear and Robson (1994) showed that the 
spectral distribution of energy in the sub-mm range of M82 
is well fitted by a grey body law with a temperature of 48 
K and an emissivity index of 1.3. 

The dust emission is, 

F v = M dust k d est B(v rest ,T)^-^- (9) 

where z is the redshift of the source, k d est is the rest- 
frequency mass absorption coefficient (e.g. Rowan-Robinson 
1986), B(u rest , T) is the rest-frequency value of the Planck 
function for dust at temperature T, D L is the luminosity dis- 
tance and (Mdust) is the total mass of dust (Hughes, Dunlop 
and Rawlings 1997). For the rest-frequency mass absorption 
coefficient (k d est ) we adopted the value of 0.15 m kg -1 as 
suggested by Hughes et al. (1997). The dependence of k d 
with wavelength is given by k d oc A _/3 , where (3 is the emis- 
sivity index of the dust grains. 

Observations of galaxies with different sources of heat- 
ing indicate an average dust grain temperature of about 50 
K (e.g. Chini et al. 1989, Chini, Kreysa and Biermann 1989, 
Hughes et al. 1993, Hughes et al. 1994). In what follows we 
assume a dust temperature of 50 K and an emissivity index 
of 1.5 given by Hughes et al. (1997). The dust mass of M82 is 
about 2 x 10 6 Mq (Hughes et al. 1994). We also compare the 
flux due to dust in galaxies hundred times brighter than M82 
(Arp 220-like luminosity, Figure 8) to the flux due to the SZ 
effect (Figure 7). We must keep in mind that Arp 220 is an 
extreme case of an ultraluminous infrared galaxy (ULIRG) 
with an infrared luminosity of about ~ lO 12 L which corre- 
sponds to 98% of the total luminosity. Combining the ob- 
served flux at 850^m (F(850^m)=792±26 mjy) given by 
Lisenfeld, Isaak and Hills (2000) with Equation 9 we esti- 
mate for Arp 220 a dust mass of about 1O 8 M0. A value of 
2xl0 8 M Q for the dust mass is given by Dunne et al. (2000) 
assuming a dust temperature of 42.2 K and an emissivity 
index of 1.2. 

Galaxies as M82 have an expected flux at 1 mm of about 
20 /iJy for redshifts between 1 and 10 and about 5 /xJy 
for observations at 2 mm (100 times higher for the case of 
an Arp 220-like galaxy). Notice that these values are sim- 
ilar to the maximum values of the fluxes due to SZ emis- 
sion for massive galaxies with baryonic masses higher than 
1O 11 M0 (Figure 7). The estimate of total fluxes due to dust, 
does not take into account the relative size of the galaxy 
since the SZ flux plotted in Figure 7 corresponds only to 
the flux estimated inside a beam size of 5". At high redshift 
(z >3) the whole galaxy will fill the 5" beam, but prob- 
ably will be resolved by interferometric observations. The 
other fact is that the dust content of M82 docs not change 
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Table 1. Main characteristics of the most sensitive telescopes at the mm range. The upper value of 10 km for ALMA correspond to the 
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maximum separation between two antennas. The mapping speed (*) is measured in aremmu e ^ 
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Figure 8. The solid lines arc the expected flux for M82 at 2 
mm (lower solid line) and at 1 mm (upper solid line) assuming a 
dust temperature of 50 K and f3 =1.5. The dashed lines are the 
expected flux for a galaxy 100 times more luminous than M82. 
The lower dashed line corresponds to fluxes at 2 mm and the 
upper dashed line, to fluxes at 1 mm. 
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Figure 9. Difference between the fluxes at 1 mm and at 2 mm. 
The different lines correspond to different galaxy masses; from 
the bottom to the top the baryonic masses are 2X10 11 5xlO n , 
10 12 , 2xl0 12 and 5x 10 12 M Q . Only the thermal SZ was included 
in this figure. The horizontal dotted lines are the estimated rms 
assuming the sensitivity and integration time shown in the upper 
left corner. 



with time, therefore the fluxes in Figure 8 can be considered 
upper limits for galaxies with lower metal content or less 
evolved than M82. 



5 DIFFERENTIAL MAPPING AND THE MM 
COLOUR-COLOUR DIAGRAM 

The particular spectral signature of the SZ effect, i.e. a pos- 
itive maximum at about 800 /im and a negative minimum 
at around 2 mm plus its non-dependence on redshift can 
be used to maximize its detection against the dust emission 
from the same source. 

Using the fact that at 2 mm the SZ has a negative 
value, the difference between the flux at 1 mm and the flux 
at 2 mm (Figure 9) is almost two times the signal given by 
the individual values (Figure 7). Notice that subtracting the 
signal at 1 mm from the signal at 2 mm will also cause a 
decrease in the flux due to sources dominated by dust emis- 
sion therefore producing a map biased to SZ sources. Source 
confusion due to dust emission should also be reduced. This 
differential mapping method will also help in removing the 
sky signal as at these low resolutions the atmosphere emis- 
sion at 1 mm and at 2 mm are correlated. In a single obser- 
vation with a wide area camera, the sky signal is removed 
by using the fact that some of the pixels are viewing the 



blank sky. Subtracting the time average level of these blank 
pixels can increase the signal to noise by a factor of 3. This 
is possible because the source structure is constant in time 
while the sky vary over the array on time scales of several 
seconds (Holland et al. 1998). The use of simultaneous mea- 
surements at different wavelengths allows to reduce the sky 
noise due to a better sampling of the fast sky fluctuations. 
Also it is possible to separate different components of the 
image (i.e. detector noise, sky, astrophysical signal) without 
assuming any atmospheric spectrum by solving a set of cou- 
pling equations and obtaining a model of the sky brightness 
distribution (Mauskopf et al. 2003). 

To distinguish between galaxies which are dominated 
by the dust emission and galaxies which are dominated by 
the SZ effect, we have explored the use of mm colour-colour 
diagnostic diagrams. 

Let us consider the emission of an object composed by 
the thermal SZ effect with a fixed comptonization parameter 
of 10 -4 plus a grey body. The amplitude of the grey body 
emission is given by changing R which is the ratio between 
the maximum of the grey body emission and the maximum 
of the SZ emission (Figure 10). 

Figure 11 represents the plane 1(2 mm)/I(l mm) vs. 1(3 
mm)/I(l mm) given by the ratios between the intensities 
at the given wavelengths. The intensity of the SZ effect is 
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Figure 10. The dash-dotted line shows the spectrum of a grey 
body (/3 = 2) at 50K. The grey body emission has been matched 
to 5 times the maximum of the SZ emission. The solid line shows 
the SZ effect with y = 10~ 4 added over the grey dust. We assumed 
a beam of 5". R. is the ratio between the maximum of the dust 
emission and the maximum of the SZ effect as defined in the text. 

fixed by adopting a y parameter of 1CF 4 . The dust emission 
is given by a grey body with temperature of 50 K and an 
emissivity of /3=1.5. The amplitude of the dust emission 
is calculated by varying R (R=0.1,l.,5.,10.,20.). The error 
bars were estimated by fixing the values of R and redshift 
but allowing the temperature to vary between 35 and 65K 
and the emissivity parameter between 1 and 2. 

Figure 11 allows us to examine the behaviour of the mm 
colours for different redshift ranges: 

• Low redshift galaxies {z 1). In this case the peak of 
the grey body is far away from the mm range where the sig- 
nal due to the SZ is dominant (between 1 mm and 3 mm). 
For R=0.1 we obtain almost the (negative) pure SZ colour 
(1(2 mm)/I(l mm)=-0.72, 1(3 mm)/I(l mm)=-0.68). For 
higher values of R the Rayleigh- Jeans tail of the grey body 
radiation starts to contaminate the SZ spectra producing a 
shift towards the right top of the 1(2 mm)/I(l mm) vs. 1(3 
mm)/I(l mm) plane, but, for these redshifts, the displace- 
ment is modest. 

• Intermediate redshift galaxies (2 ^ z Si 5). This case is 
similar to the previous one but the grey emission on top of 
the SZ is more important, producing a larger displacement 
of the points to the right top corner of the diagram, while 
the mm colours remain negative. 

• High redshift galaxies (z > 5). In this case the R-J part 
of the grey body emission has covered the SZ effect emission 
at 2 mm, producing in some cases positive 1(2 mm)/I(l mm) 
and 1(3 mm)/I(l mm) ratios. 

Figure 11 indicates that all points in the bottom left 
part of the diagram have a dominant SZ component. This 
diagnostic seems valid up to redshift of several. The accuracy 
or reliability of the diagnostic will clearly depend on the size 
of the error bars that in turn depend on the amplitude of 
the signal and sensitivity of the telescope. All in all it seems 
that the new facilities would permit to start the exploration 
of at least the largest star-forming systems at z < 12 even 
if they show moderate dust emission. Inside the estimated 
error bars, sources dominated by the SZ effect fall clearly 



in the distinct quarter of the diagramme where both the 
2 to 1 and 3 to 1 mm intensity ratios are negative. This 
result is very robust in the sense that sources dominated 
by the SZ are going to be separated from the dust emission 
independently of redshift. The confusion is higher for high 
redshift galaxies with the dust emission comparable to the 
SZ emission. 

The diagnostic diagram presented in Figure 11 would be 
partially covered by future deep surveys at mm wavelengths, 
showing high redshift sources dominated either by dust or 
by SZ emission. Because of the expected delay between the 
burst event and the production of metals and dust, our di- 
agnostic diagramme can potentially be used to constrain the 
sources and time scales for dust and metal production in the 
early universe (Morgan & Edmunds 2003). 

To further advance in the study of the evolution of the 
SZ and dust emission in a self-consistent manner, to check 
in more detail the possibility of observing the SZ through in- 
dividual galaxies with the next generation of mm telescopes 
and to estimate the contribution of SZ sources to the back- 
ground radiation, we will present in a forthcoming paper 
a set of simulated mm maps including sources dominated 
by SZ and those dominated by dust emission (with different 
clustering properties) plus the comparison of our predictions 
on the SZ background with those made by Aghanim, Bal- 
land & Silk (2000) based on the SZ due to black holes-seeded 
proto-galaxies. 



6 CONCLUSIONS 

In this paper we have explored the possibility that the cen- 
tral comptonization y parameter in young spheroidal galax- 
ies or galaxies undergoing a massive starforming event may 
reach values comparable to those observed in galaxy clus- 
ters. 

However due to the poor angular resolution-sensitivity 
combination of present day mm facilities, the detection of 
the SZ effect in individual galaxies may only be possible 
with the next generation of high angular resolution - high 
sensitivity mm telescopes (e.g. GTM/LMT or ALMA) and 
then probably only in the most massive systems capable of 
producing a 800 fxrn flux of more than a few /xJy for periods 
of time of several hundred megayears. 

We have included instrumental noise in our calculations, 
but neglected the noise contribution caused by high red- 
shift millimeter background point sources (e.g. Blain 1998, 
Hughes et al. 1998). The confusion due to unresolved sources 
depends on the actual details (i.e. resolution) of the observa- 
tion, the (unknown) emission radial profile of these sources 
and the extrapolation of the number counts to fainter fluxes. 
Blain et al. (2000) give an estimation of the confusion limit 
of about 10 fj,3y for the case of observations at 1 mm and a 
beam size of 5". This value is comparable to the expected SZ 
flux of galaxies with baryonic masses of 2xlO 12 M0 and an 
age of ~ 3 x 10 8 years (Figure 7). Confusion will be a major 
problem when combining independent observations at dif- 
ferent frequencies but its effect will be greatly reduced using 
instruments that simultaneously observe at three selected 
frequencies. 

A major difficulty in the detection of the SZ effect in an 
individual galaxy will be its confusion with the dust emis- 
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Figure 11. The different panels represent a set of mm colour-colour diagram for sources with both the emission of the SZ effect and 
the emission from the dust. The different panels correspond to galaxies at different redshifts: z=0.1 (top left), 2=3 (top right), z=8 
(bottom left) and 2=15 (bottom right). For each individual panel, the different points correspond to the adopted values of the ratio R 
=0.1,1. ,5. ,10. ,20. The amplitude of the bars are calculated by varying the grey body temperature from 35 to 65K and the emissivity 
parameter from 1 to 2. The lines correspond to T=50K grey body models with emissivity parameters 1 (solid line), 1.5 (dotted line) 
and 2 (dashed line). The values corresponding to z=12 and 2=20 are indicated. The cross in the bottom right corner corresponds to 
observational errors for 3cr detection in all three bands. 



sion originated in the same star formation event. To solve the 
dust contamination problem we suggest to exploit the par- 
ticular spectral signature of the SZ spectrum and maximize 
its segregation from dust emission by using simultaneous 
multifrequency observations or differential detection at 1, 2 
and 3 mm. In the diagnostic diagram presented in Figure 1 1 
there is a wide region centred at 1(2 mm)/I(l mm)=-0.3 
and 1(3 mm)/I(l mm)=-0.3 where the thermal SZ emission 
can be clearly separated from the dust emission over a wide 
range of redshifts. 

Due to its particular spectrum the detection of the SZ 
signature in individual galaxies and the discrimination be- 
tween sources dominated by the SZ effect from those sources 
dominated by dust emission would be enhanced by the use 
of differential mapping systems that simultaneously observe 
at three selected frequencies (Section 3). Additional advan- 
tages of such systems will be the lowering of the minimum 
detected flux due to a better sampling of the time corre- 
lated atmospheric fluctuations (i.e. reducing the sky noise) 
and the relative absence of source confusion generated by 
dust emitting sources. 

New generation of mm cameras such as SCUBA-2 t 

t http: \\www.roc.ac.uk\atc\projects\scuba_two 



or SPEED i - which will have the capability of observing 
simultaneously in more than a band - would be ideal to 
detect compact and faint SZ sources. 
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